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Abstract 

On the basis of the density functional theory, we give a clear definition of 
the free energy landscape. To show the usefulness of the definition, we construct 
the free energy landscape for rearrangement of atoms in an FCC crystal of hard 
spheres. In this description, the cooperatively rearranging region (CRR) is clealy 
related to the hard spheres involved in the saddle between two adjacent basins. A 
new concept of the simultaneously rearranging region (SRR) emerges naturally as 
spheres defined by the difference between two adjacent basins. We show that the 
SRR and the CRR can be determined explicitly from the free energy landscape. 
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1 Introduction 



The thermodynamics and the dynamics of glass forming substances exhibit anomalous 
behaviors near the glass transition temperature, T g . As the sample which was prepared 
by rapid cooling below T g is heated, the specific heat exhibits an abrupt increment 
at T g [1]. It is also known that supercooled liquids near T g show several relaxation 
processes and that the relaxation time of the slowest process tends to diverge at a 
temperature below T g [2] . It is, therefore, desirable to construct a theory for the glass 
transition which gives a unified explanation for these dynamic and thermodynamic 
anomalies. 

The mode coupling theory (MCT) [3] for the glass transition describes the liq- 
uid dynamics based on the generalized Langevin equation and predicts an ergodic- 
to-nonergodic transition at a critical temperature, T c due to the nonlinear coupling 
between various modes. However, it is now believed that T c is much higher than T g . In 
addition, MCT cannot treat the thermodynamic singularities. Mezard and Parisi esti- 
mated thermodynamic properties of glassy state in quasi-equilibrium with the replica 
method [4]. In this approach, the Kauzmann temperature, Tk, becomes an ideal glass 
transition temperature. However, Tk is much lower than T g and the replica method 
cannot treat the dynamics either. Both theories thus can not give a unified explanation 
for the thermodynamic and dynamic singularities of supercooled liquids near T g . 

We have proposed the free energy landscape picture by which a unified understand- 
ing is given for the singularities near T g [5]. The free energy landscape is defined 
in the configurational space as a function of fictitious localized atomic positions and 
each point in the free energy landscape is determined by coarse-graining of the mi- 
croscopic motion in the configurational space. The free energy landscape is expected 
to have many basins and the saddle points at low temperatures. The separation of 
the dynamics can be understood by the separation of the dynamics among basins and 
the dynamics within a basin. The thermodynamic and dynamic singularities near T g 
can also be explained by the free energy landscape picture. The trapping diffusion 
model [6] and its generalization to dynamics in a landscape [5] provide a unified under- 
standing of the dynamical singularities near the glass transition. This model showed 
that the dynamical transition occurs at the temperature where the mean waiting time 
for escaping from a basin diverges. Tao et al. [7] proposed a model based on the free 
energy landscape picture to explain the specific heat anomaly at the glass transition. 
They showed that the anomaly of the specific heat occurs when the mean waiting time 
exceeds the observation time at low temperatures. 

It is now believed that the glass transition can be phenomenologically well explained 
by the free energy landscape picture. It is thus an important question how one can 
explicitly define and construct the landscape on the basis of the microscopic Hamilto- 
nian. Because of this lack of an explicit theory, the analysis of the dynamics in the 
landscape has been detered so far. 

The aim of the present paper is to define explicitly the free energy landscape using 
the microscopic Hamiltonian. To this end, we exploit the density functional theory 
(DFT) [8-11]. 

One of the key concepts to understand the glass transition singularity is the cooper- 
atively rearranging region (CRR) proposed by Adam and Gibbs [12]. Adam and Gibbs 
succeeded in explaining the non-Arrenius behavior in the temperature dependence of 
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the viscosity for fragile glass forming liquids. However there is no clear common un- 
derstanding of the physical meaning for the CRR. From the present definition of the 
free energy landscape, the clear meaning of the CRR emerges. In addition to the CRR, 
the present construction of the landscape leads to a new concept of rearranging region, 
namely simultaneously rearranging region (SRR) which is defined by the difference 
between the atomic configuration in two adjacent basins. 

This paper is organized as follows. We give a definition of the free energy landscape 
by the DFT in Sec 2. To confirm the usefulness of this definition, we construct the free 
energy landscape for rearrangement of atoms in an FCC crystal of hard spheres in Sec 
3. A relation between the free energy landscape and the CRR is given in Sec 4. We 
also discuss the SRR in Sec. 4. Results are summarized in Sec 5. 

2 Definition of the free energy landscape 

We define the free energy by the partition function which is obtained by the partial sum 
of the fast motion in the configurational space. The free energy becomes a function 
of the fictitious atomic positions and we call this free energy as the landscape. We 
exploit the density functional theory (DFT) [8-11] to calculate the free energy. In 
the DFT the grand potential, f2[p(r)], is expressed as a functional of the density field, 
p(r), and the equilibrium density field, p e(? (r), is determined as a minimum of the 
grand potential. If the density field is given as a function of a particle configuration, 
{Rj} = {Ri, R2, • • • , Rjv} where R; is the position of the i-th particle, the density 
functional f2[p(r)] becomes a function of {Rj} which can be considerd as the free 
energy landscape defined above. We employ a sum of Gaussians as the density field: 



Here a and Rj are the degree of the spread of the density distribution and the position 
of the i-th particle, respectively. This density field has been used for the investigations 
such as the liquid-solid transition [9,10,13,14] and the glass transition [15-19]. The 
density field (JIJ means that the distribution of the vibrational motion around {Rj} is 
approximated by Gaussian functions. The vibrational motion within |r — Rj| < \j \fa 
is coarse-grained by a. By using eq. (JTJ), f2[p(r)] becomes a function of a and {Rj}: 



We use the grand potential (J2J) as a definition of the free energy landscape. Our 
definition of the free energy landscape is different from the definition by Dasgupta and 
Vails [20]. They defined the free energy landscape as a functional of the density field, 
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(1) 



nD>(r)] = fi(a,{Ri}). 



(2) 



P(r). 
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3 Calculation of the free energy landscape 
3.1 Model 

We use the approximation proposed by Ramakrishnan and Yussoff to calculate f2[p(r)] 
[21]. The result of the approximation is as follows: 

(3An[ P (r)} = pn[ P ( r )]-pn[ Pl ] 

= [ drp(r)log[^l - / dv{p{v)- Pl ) 
Jv 1 pi J Jv 

j dv x J dv 2 c{\r l - r 2 |)(p(r 1 ) - pi){p{r 2 ) - pi). (3) 

Here (3 = (ksT)~ l where ks is the Boltzmann constant and T is the temperature, 
and pi denotes the density of the uniform liquid. The coefficient c(|r|) is the direct 
correlation function of the uniform liquid at p\. 
We substitute eq. (JTJ) into eq. (jUJ): 

/3Ml(a,{Ri}) = F, d (a,{R i }) + F -^/o(«)-iEE / K R u) (4) 



where 



F ld (a,{R l }) = (-) i Yl I rfrexph^r-R^logiir^^exph^r-R,) 2 ]! 

7T j Jv [P l 71 j J 

(5) 

F = N l (p 8 -%) I dvc{r) (6) 



2 J Jv 



I (a) = (-) / dr x [ dr 2 c(|ri - r 2 |) exp[-ar 2 ] exp[-ar 2 ] (7) 
v 7r y Jv Jv 

J(a,R) = f-) 3 / dri [ dr 2 c(\r 1 - r 2 |) exp[-ar 2 ] exp[-a(r 2 + R) 2 ]. (8) 

V 7T 7 JV JV 

The average density p s is defined by 



p s = I / drp(r), (9) 



and JV a = p,V, iV, = p z F. 

In the present paper, we treat a system of the hard spheres. We use Percus-Yevick 
approximation for the direct correlation function of the hard sphere [8]: 

= + % + V (r<a) 

w [ (r > a). y ' 

where a is the diameter of the hard sphere. The coefficients B , Bi, and B 2 in eq. (|1U|) 
are given by 

» {l + 271)2 (U) 

(l + rp/2) 2 

Bi = 6 ^7r^r (12) 

B, = -*£±*£ , (13) 
2 (1 -rp) 4 ' v ; 
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and rj = npi/6. By using eq. (jTUjl . F , Jo (at), and I(a, R) can be calculated analytically 
[14,22]. 

We use aa 2 = 478 correponding to the liquid-solid transition point [13]. For large 
a, Fid(a, {R*i}) reduces to 



F„(a,{R,}) - iV.ftl/drexph^logr^t^lp)^ 

_l og Q_lo g ( A )__] +p; . ( M ) 



Though Jones and Mohanty proposed that Jo (at) and J(a, R) also reduce to more 
simple form [14] for large a, we do not use their expression, since it can be used only 
for large a, and is not applicable for the entire region of a. We use the analytical form 
which can be used for all a [22]. 

3.2 Forced relaxation in an FCC lattice 

In order to show the usefulness and the validity of the present definition of the free 
energy landscape, we calculate the free energy landscape for a forced structural change 
in an FCC crystal of hard spheres. Since the grand potential (J2J) has too many parame- 
ters, a and {Hi}, we calculate the free energy landscape for a process in which we force 
to move a specific group of particles and to relax other particles to minimize Q(a, {Rj}) 
in the {R-i} space. We calculated the free energy landscape when we force to rotate 
two adjacent particles with their distance fixed (Figure 1). The axis of rotation is the 
center of the two particles and two particles are rotated around an axis in < 100 > 
direction. The arrangement of other particles is changed to minimize Q(a, {Ri}) in 
the {Rj} space with the positions of the two particles fixed. We exploit the steepest 
decent method to find the minimum. 

We imposed periodic boundary conditions in all directions. The number of the 
particles in the system is 108. We set p s a 3 = 1.0 and pi<7 3 = 0.967 which correspond to 
the metastable state to the liquid state [13]. Figure 2 shows the free energy landscape 
as a function of rotation angle 9. The free energy landscape has a saddle point at 
9 = 7i / 2 and two minima at 9 = and 9 = tt. Two minima at 9 = and 9 = tt 
correspond to the FCC configuration (Figures 3 (a) ~ (e)). The height of the saddle 
point is about 57ksT '. 



4 CRR and SRR in the free energy landscape 
4.1 Definition of CRR and SRR 

We define the simultaneously rearranging region (SRR) as the difference in the config- 
uration at the minimum between the two adjacent basins which are connected directly 
via a saddle point. Apparently, the number of the particles in SRR is 2 in the system 
studied in Sec 3. 

The CRR corresponds to the particles which are moved at the saddle point from 
the configuration at the minimum of the original basin. Our definition of the CRR is 
thought to be closely related to the idea of the CRR by Adam and Gibbs [12], since the 
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particles which are moved at the saddle point are thought to give the configurational 
entropy for the rearrangement proposed by Adam and Gibbs [12]. 

4.2 Calculation of the CRR 

Since we use aa 1 = 478, the degree of the coarse-graining is about 0.05cx. Hence we 
assumed that the particle in the CRR for the system calculated in Sec 3 is the particle 
displacing more than 0.05<r at 9 = 0.5n from the position at 9 = 0. 

We found that the number of the particles in the CRR is 46. The particles which are 
shown by the bold circles at 9 = ir/2 in Figure 3 are those which satisfy the criterion. 

5 Summary 

In the present paper, we have given a clear definition of the free energy landscape using 
the density functional theory. As a test of the definition of the free energy landscape, 
we have calculated the free energy landscape for rearrangement of atoms in an FCC 
crystal of hard spheres. 

We have clarified the concept of the cooperatively rearranging region (CRR) by 
giving a relation between the free energy landscape and the CRR. We also have intro- 
duced a new region, simultaneously rearranging region (SRR) which is different from 
the CRR. We have calculated the number of the particles in the SRR and the CRR 
from the calculated free energy landscape. 

The present frame work can be easily generalized for glass-forming substances. 
Hence we can calculate the number of the particles in the CRR and the height of the 
saddle point near the glass transition from the free energy landscape. Details will be 
studied in a forthcoming paper. 
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Figure Captions 



Fig 1 A schematic representation of the model under consideration. Two particles are 
rotated with their distance fixed. The arrangement of other particles are changed 
to minimize Q(a, {Rj}) in the {Rj} space. 

Fig 2 The free energy landscape when we force to rotate the two particles. The 
horizontal axis is rotation angle 9. The vertical axis is the difference between 
n(a, {Rj}) of a certain rotation angle and Q(a, {Rj}) of FCC. The number of the 
particles in the system is 108. We set p s cr 3 = 1.0, pier 3 = 0.967, and aa 2 = 478. 
Although we write the horizontal axis by one dimension, the actual free energy 
landscape has 108 x 3 dimensions. 

Fig 3 The configurations at the plane which two particles rotate. Each figure corre- 
sponds to (a) 9 = 0, (b) 9 = 0.3tt, (c) 9 = 0.5tt, (d) 9 = 0.7tt, and (e) 9 = ir. 
Two particles with the bold circle are the particles rotating. The particles with 
the bold circle at 9 = 0.5n are the particles in CRR. 
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Fig 2 
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Fig 3 
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